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^ ! Abstract 

We continue to study the effect of uneven voltage distribution on two close cylindrical conductors 
with parallel axes started in our papers fl and Q], now to find the electrostatic torques. We 

o', 

Q" 1 calculate the electrostatic potential and energy to lowest order in the gap to cylinder radius ratio for 

an arbitrary relative rotation of the cylinders about their symmetry axis. By energy conservation, 
the axial torque, independent of the uniform voltage difference, is found as a derivative of the 
> i energy in the rotation angle. We also derive both the axial and slanting torques by the surface 



integration method: the torque vector is the integral over the cylinder surface of the cross product 
of the electrostatic force on a surface element and its position vector. The slanting torque consists 
of two parts: one coming from the interaction between the patch and the uniform voltages, and the 
other due to the patch interaction. General properties of the torques are described. A convenient 
model of a localized patch suggested in Q] is used to calculate the torques explicitly in terms of 



X 

elementary functions. Based on this, we analyze in detail patch interaction for one pair of patches, 
namely, the torque dependence on the patch parameters (width and strength) and their mutual 
positions. The effect of the axial torque is then studied for the experimental conditions of the 
STEP mission. 
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I. INTRODUCTION 



The actual distribution of charges on a metal does not guarantee its surface to be an 
equipotential because of the impurities and microcrystal structure of the material. This 



phenomenon, known as patch effect (PE; for its experimantal study see paper is re- 
sponsible for the mutual force and torque between two metallic surfaces at finite distances. 
The effect is the larger, the closer the surfaces, as first confirmed by the calculation of the 



patch effect force for two parallel conducting planes 14 1. 

n n 

PE in a cylindrical geometry was studied in the first two parts of our paper, 

\M and M 

(henceforth referred to as CPEI and CPEII, for '"Cylindrical Patch Effect"'), where the PE 
energy and force have been examined. Here we calculate the torque due to PE between two 
coaxial cylinders. This calculation completes the study of the electrostatic interaction for a 
cylindrical capacitor; our analysis largely benefits from the results found in CPEI and II. 

PE is important for any precision measurement if its set-up includes conducting surfaces 
in a closed proximity to each other. For the STEP experiment [5 8 ] , where the differential 
axial motion of cylindrical test masses (TM) will be used to test the universality of free fall 
to an unprecedented accuracy of about 1 part in 10 18 , the axial torque is of a particular 
interest (see section IVI]) . 

We determine PE torques between the two boundary surfaces of an infinitely long cylin- 
drical capacitor. By energy conservation the rotation by an angle 7 about the direction 
7 of one of the conductors relative to the other causes an electrostatic torque in the same 



direction which is given by the formula (see, for instance, {9]) 

where Wi^y) is the electrostatic energy as a function of the rotation angle. However, due 
to the specifics of cylindrical geometry, we can properly imagine a rotation only about the 
symmetry axis of the two infinite cylinders. Tilting, say, the inner one about any other 
direction leads to the intersection of the cylinders at some finite distance, i. e., to the 
breaking of the problem geometry. For this reason, we employ also a different method of 
the torque calculation. The force, dF, due to the electrical field E acting on a small area 
dA of a conductor with the charge density a is given by 

dF = aEdA , (2) 



(see our comment [1CJ). The resulting element of the torque about the origin at distance r 
from dA is then 

df = fxdF , (3) 

which expression, integrated all over the body surface, gives the general expression of the 
torque acting on the conductor. 

The energy, the field and the surface charge density which are needed in the formulas ([1]) 
and ([3]) are expressed through the electrostatic potential in the gap. For typical experimental 
conditions, such as in the STEP configuration {(3, 7], the gap, d = b — a, is much smaller 
than either of the cylinder radii, a < b. This, first, justifies the model of infinite cylinders, 
especially if the patches are predominantly far from the real cylinder edges, and, second, it 
allows for a significant simplification of results to lowest order in dja. 

In the next section we solve the boundary value problem (BVP) for the potential in the 
gap with general voltage distributions on the cylinder surfaces. Based on this, we find the 
energy in section [TTTl and then the longitudinal PE torque by the formula ([1]). In section HVl 
we derive all the three components of the PE torque by the surface integration method using 
formula ([3]) (the two expressions for the axial torque agree precisely). In section |V] a model 
of the localized patch potential introduced in CPEII is described, ending with closed-form 
expressions for the torques. The latter are then calculated and analyzed in the case when 
a single patch is present at each of the cylinders. In section [VI] we apply our results to the 
STEP experiment set-up, coming up with a clear picture of the test mass rotation motion. 
The details of calculations are given in two appendices. 



II. ELECTROSTATIC POTENTIAL 



We employ both Cartesian and cylindrical coordinates in two frames related to the inner 
and outer cylinders as shown in fig. [U In the outer frame an arbitrary point is labelled by 
the vector radius r , Cartesian coordinates {x , y , z}, or cylindrical coordinates {p, if ', z}; 
in the inner frame the corresponding quantities are f, {x, y, z}, {p, if, z}. The origins of 
the frames coincide so that the primed and unprimed coordinates are related simply by a 
rotation, by some angle 7, about the z axis: 

x = a; cos 7 + ysiaj, y = — a; sin 7 + ycosj, z —z, (4) 
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or, in cylindrical coordinates, p' = p, = if — 7, z' = z. The surfaces of the inner and 
outer cylinders are thus described by the equations p = a and p = b, respectively, and are 
assumed to carry arbitrary distributions of electrostatic voltage. Hence the electrostatic 
potential, $, satisfying the Laplace equation in the gap between the cylinders, 



A$ = 0, 



p > a, p < b, < if < 2n, \z\ < 00 



satisfies also the boundary conditions of the first kind at the cylinder surfaces: 



G((p,z) ; 



p=a 



p=b 



V- + H(<p', z) = V~ + H(<p - 7, z) 



(5) 



(6) 



Here V~ = const is the uniform potential difference, so all voltages are counted from the 
uniform voltage of the inner cylinder taken as zero. The non-uniform potential distributions, 
i.e., the patch voltages, are described by arbitrary smooth enough functions G(tp,z) and 
H(tp,z). Same as in CPEI, II we assume these functions squarely integrable; wherever 
proper, we will also assume them, as done in conditions (10), (All) (A12) and (C5) of 
CPEI. 

For any squarely integrable function u((p, z) we have its Fourier expansion with the Fourier 
coefficient u n {k): 

2-n 00 

u(<p,z) = — j_ dk u n (k)e likz+n ^ , u n (k) = — I I d^dzu(ip,z)e-^ kz+n ^ . (7) 



2tt 



-00 



For any two such functions u((p, z) and v(ip, z) the useful Parceval identity holds, 



2n 00 



[U, V = 



J d(pdzu((f,z)v*((f,z) = J dk ^ u n (k)v*(k) ; 



(8) 



-00 



here and elsewhere the star denotes complex conjugation. 

According to the boundary condition ([6]), we split the potential in two parts due to, 
respectively, the uniform boundary voltages and patches: 



$(f) = $ u (f) + $ p (r) 



p=a 



V 



i{kz+nip) 



I poo 00 

= G(<p,z) = — / dk £ G n (k)e 

a 111 J-00 n= _ oc 

1 roo 00 

= H(<p- 1 ,z) = — dk H n (k)e- m ^e 

, Z7T J-00 „_ „ 



(9) 
(10) 

(11) 



iwy i(kz+mp) . 



Function $ M is the classical solution for a cylindrical capacitor [9]; to 1. o. in d/a it is 
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$«(r) = (a/d)V- ln(p/a) 



(12) 



Function <3> p is gotten by the standard separation of variables in cylindrical coordinates [see 



c.f. 



ll| . Chs. 5, 6]. Its representation satisfying formally the Laplace equation is: 



1 PCX) °° 

® P (r) = 7T dk E \Mk) 4(M + B n (k) K n (kp)] e 

Z7T J-oo „ — 



i(kz+nip) 



(13) 



where I n (0> -^n(C) are the modified Bessel functions of the 1st and 2nd kind, respectively 
[the Macdonald's function K n (£) definition for the negative values of its argument is taken 
by the parity of so > K n {kp) stands for (sign k) n K n {\k\p)}. The unknown A n (k), B n {k) 

are found form the linear system implied by the boundary conditions (jlll) : 



A n (k) I n (ka) + B n {k)K n {ka) = G n {k) ; 
A n (k)I n (kb) + B n (k)K n (kb) = H n (k)e~ in \ 



(14) 



n = 0,±l,±2, ... . 



It is the same system that has been effectively solved, to lowest order in d/a, in CPEI 
[Appendix A, coefficients A^(k) and B®(k)], with the exception of e~ ini in the r.h.s. instead 
of e lhz ° . So, by replacing H n (k)e~ tkz ° with H n (k)e~ tn ' 1 in the answer (A13), CPEI, we get: 



A n (k) = -- d {K n {kb) [G n {k) - H n {k)e~^] } , (15) 
B n (k) = ~ {l n (ka) [G n (k) - H n (k)e-^} } . 

Thus the electrostatic potential ([13]) . to 1. o. in d/a, is: 

$ p (r) = -- / dk \Gn(k) - H n {k)e-^\ Q n (kp)e< kz+n ^ ; (16) 

fi„,(A;p) = K n (kb)I n (kp) - I n (ka)K n (kp) . (17) 
Formulas f|T2|) . f|T6|) . and ( TT7|) allow us to calculate both the electric energy and field. 



III. AXIAL TORQUE BY THE ENERGY METHOD 

The uniform potential ( 1T2|) does not depend on 7, so the variation of the electrostatic 
energy due to the rotation comes only from the patch potential ( JT6l) . same as it happens with 
the axial PE force (CPEII, section III). The axial torque, T z , is thus given by the formula 
JT]) where W{^f) is replaced with W p {^) that we calculate below. 
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A. Electrostatic Energy 



Denote T>oo the infinite domain between the two cylinders of our capacitor. The patch 
energy stored there, finite due to the locality of patch distributions, is: 



W p = - f (V$ p ) 2 dV = 

!roo p2ir 
/ dz bd(p H(ip — 7, z) 
J-oo Jo 



£o 
2 



dp 



P =b 



oo r2TT 

dz / ad(p G((p, z)— — 
-oo Jo op 



(18) 



p=a 



Here we used boundary conditions f ill I) and the fact that the potential is harmonic in the 
domain l^, see section III in CPEI for details. The double integrals above are calculated 
via Fourier coefficients of the potential and its derivative in p by the Parceval identity (JHJ). 
The Fourier coefficients of the derivatives are found from the formula f|T6|) ; the calculation 
goes the same way as in CPEI, Appendix C, and results in 



dp 



p=a,b 



1 roo °° r 

■dLoo^ £ [G n (k) - H n (k)e 



■i n-y 



(19) 



To 1. o. in d/a, this expression holds at both the inner and the outer boundary. Using the 
Fourier coefficients, G n (k) and H n (k), of the boundary functions, we write formula (I18p as: 

(20) 



W? = ^ Loo dk ^ \ Gn{k) ~ Hn{k)t 



-in-y 



The only part of this that depends on 7, and thus contributes to the axial torque, is: 

W p { 7 ) = — °- / dk ^\G n (k)H* n (k)e m ~< 



(21) 



B. Axial Torque 

Using (12"T|) . we calculate the axial torque by the formula (jTJ): 

dW p (a) 



T, 



' d 



dk 



n S 



G n (k)H* n (k)e 



m-y 



(22) 



This representation is valid, to lowest order in d/a, for an arbitrary rotation 7. The torque 
does not vanish only if patches are at both boundaries [G n , H n ^0]. A non-zero torque 
is generally found for 7 = 0, unless G n {k) = XH n (k), X real, i. e., the patch distributions 
at both cylinders are the same up to scaling. The expression of the axial torque perfectly 
matches that of the axial force in the symmetric configuration [CPEII, formula (21)]; the 
only difference is 7 instead of z°, and the factor n in place of k. 



IV. ALL TORQUES BY THE SURFACE INTEGRATION METHOD 



A. General Formulas for Electrostatic Torques 



According to the formula ([3]), the patch effect torque on the outer cylinder is: 



f = b f 277 dp f°° dzaffx E) 

J J — oo 



(23) 



r2-K 

e a I dip 
o 



dz 



P =b 
<9$ 



r x 



19$ 



dp p dp 



dz 



dp 



P =bi 



\ a , 



Here we expressed the electrical field through the potential, and the charge density as the 
product of eo and the normal component of the field, by the Gauss law. We also set 
b = a + d ~ a, so the above holds to 1. o. in d/a. The Cartesian components of the torque 
are found using the well-known cylindrical unit vectors: 



6od I dip 




dz 



zd<$> 
pdp 



<9$ d& 



cosp 



p 



zd§ ( <9$ <9$\ 

pd^ s[nip+ { p dI- Z ^ 008 * 



dz 



y + 



z- smtp 



x 



(24) 



'd& 



d<$> 



dp dp 



P =b 

Because of the bilinear structure of this expression, the splitting ([9]) of the potential in the 
sum of $ u and $ p implies formally three contributions to the torque: one from the uniform 
voltages only, the other due to the interaction of patches and uniform voltages, and the 
third one from the patches only, just like we had it for the force in CPEII. However, uniform 
voltages do not give any torque, so the first contribution vanishes. Likewise, the interaction 
between the patches and uniform potential difference gives zero axial torque: the interaction 
torque is perpendicular to the symmetry axis (slanting torque): 



T 



Int 



e a- 



-\Irit 



dp 

d$ u 



P =b 



2tt 



dp 



dz —2z— — sin <p — cos p> + p— — sin cp 

oo \ dp ' p dp ' dz , 



€qCI- 



dp 



P =b 



2 71 



dp 



oo / Q$P 

dz I 2z— — cos p> 



zd<$> p 



dp 



P dp 



smp — p- 



dz 



cosp 



P =b 



P =b 



(25) 



(26) 



In contrast with that, the torque due to the patch interaction generally has all the com- 
ponents: 



rpp 

X 



r 2ir 

eoa / dp 
o 



d<5> 



p 



oo^ I dp 



z <9$ ? 



P d V 



cos p + p- 



dz 



d® p \ 
z— — I snap 



dp 



(27) 



P =b 
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Formulas (I25j) — ( 1291) provide the general representation for the torque based on the surface 
integration method. Below we use them, along with the expressions ( TT2l) and ( |T6|) of the 
potential in the gap, to find the PE torque on the cylinder. 



B. Axial Torque by the Surface Integration Method 



We compute the axial component ([29]) of the torque employing the Parceval identity (|8]): 

1 



/OO " v 
dk £ 
-OO 



2d 



G* n (k)-H*(k)e m ^) mH n (k)e- m ^+ 



(30) 



^ (G n (k) - H n {k)e-^) inH* n (k)e in i 



~ r dk ]T 3 \nG n {k)H*(k)e in ^' 

" J ~°° n=-oo 



We used Fourier coefficients of the two derivatives of the potential, the radial one f ll9p . and 
the angular one obtained by differentiating the second of boundary condition ( TTTT) : 



d<$>P 



dip 



1 /"OO °° 

= — / dk e i{kz+nv) 

p=b 27V J -°° n=-oo 



inH n (k)e~ 



-1717 



(31) 



The axial torque f[3"Tj) and (J2"2~j) by the surface integration and energy method, respectively, 
is exactly the same; this is an important cross-check of our calculations. 



C. Slanting Torques by the Surface Integration Method 



1. Uniform and patch potential interaction 



We start with calculating the torque due to the interaction of uniform and patch poten- 
tials. We first substitute the uniform field, V~/d, in the formula (j2"5|) : 

l x — — —V I dip / dz \—2z — — smip — cos </? + p— — sm tp 

d Jo j-oo y dp p dip dz J 

to get the x component. The last term vanishes after integrating in z, so: 



P =b 



T 



Int _ fO^ 

d 



rl'K rOO I 

V~ / dip dz \ -2z 

JO J-oo \ 



dp 



sm cp — 



zd$ p 
p dip 



COS ip 



p=b 



By the definition flTJ) of the Fourier transform, this double integral is equal to 2tt times 
the Fourier coefficient of the integrand at n = k = 0. The needed Fourier coefficients are 
determined in Appendix A, formulas (1A7j) and (1A8[) . so: 



T 



Int 



or 



^-(G^-H^ky 



dk 













HQ] 




k=0 





(32) 



This final compact formula for the torque has been obtained by employing the property 
G- n (—k) = H- n (—k) = H*(k) of Fourier coefficients of real functions. The estimate 

of the remainder in the expression (132]) holds for patch distributions G((p, z) and H(<p, z) 
satisfying conditions (All), (A12) and (C5), CPEI, and also such that the products zG(<p, z), 
and zH((f, z) are squarely integrable, see formula (lAlOj) in Appendix A. The validity of these 
conditions is assumed everywhere below, including the final expressions of all the torques. 

A similar calculation for the y component starts from the formula ( 126]) and uses the 
Fourier coefficients (]A6j) and (]A9j) . It results, to lowest order in d/a, in: 

() (Gx(k) - H x {k)e 



rplnt 

V 



4:71— — V ^ 

d 2 



-17 



dk 



(33) 



fc=0 



with the same remainder as in formula 



2. Patch potentials interaction 



Now we go for the expressions of the torque caused by the interaction between patches. 
The x component of this torque, TJ", given by the formula ( 12 7p . with the help of the Parceval 
identity becomes [the proper Fourier coefficients are found in (1191) . ( 1A7I) . ( 1A8j) . and ( 1A12I) ]: 

e a roa 



TP 

X 



/oo — _ 
dk £ 
-oo __ 



G* n (k) - H* n (k)e 



X 



d_ 

dk 



\ a . 



2t/^ J — oo 

(G n ^(k) - H n _ x {k)e- i{ - n -^) - (G n+1 (k) - H n+1 (k)e^ n+1 ^) 

This formula can be simplified further: integrating by part in the first of the two products and 
shifting there the index n by one, n' = n — 1, leads to the final more compact representation: 

d 



T * = w r dk £ ^1 \ G *^ - E ^y n " 



dk 



(k)-H n+1 (k)e 



-i (n+l)7 



• (34) 



The other component, T£, can be determined in a similar way starting with the expression 
( 128]) and combining it with the formulas (JT9J) , (jA"6]l . dA9]) . and (IMip . The result is: 



y 



-— r dk v 

G?^ J —oo „ ^, 



G* n {k) - H*{k)e 



d_ 

dk 



(k)-H n+1 (k)e 



-i (n+l)7 



, (35) 



entirely similar to ( )34|) . The general analysis of PE torques is completed. 
D. General properties of the PE torque 

Looking at the results of the calculation of the PE torques, one can come up with a few 
general conclusions regarding their properties, such as: 

1. Uniformly charged cylinders do not give rise to any torque. 

2. The axial torque is inversely proportional to the gap width, the transverse components 
are inversely proportional to its square. 

3. Patches need to be present on both the cylinders to generate an axial torque. 

4. A non-zero axial torque is generally found when the cylinders are not rotated against 
each other (7 = 0), unless the patch voltage distributions on both cylinders are the same up 
to a factor. 

5. Just one patch is enough to generate a slanting torque. 

6. A non-zero slanting torque is generally found when 7 = unless the patch voltage 
distributions on both cylinders are the same. 

7. The interaction between patches and uniform potentials involves only the first har- 
monics of the azimuthal angle of the patch distribution. 

So, the general formulas that we obtained enable one to make some significant conclusions 
about PE torques agreeing with the physical insights into their origin. 

V. SINGLE PATCH AT EACH OF THE ELECTRODES: A PICTURE OF PATCH 
INTERACTION 

To study the features of PE torques we analyze them in the case when only one localized 
patch is found at each of the cylinders, as it was done with PE force in CPEII. To make the 
analysis results transparent, one needs to have the torques in a simple enough closed form, 
which requires some special choice of the generic patch model, a rather delicate task. The 
localized potential distribution that satisfies this very well has been suggested and developed 
in CPEII, section IV. We repeat basic facts about this model and use it to calculate the 
torques and examine the patch interaction. 
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A. The Patch Model 



The suggested model of the patch potential is: 



V(<p -ip*,z - z„) = V((f - (p*, X,z- z„, Az) = V* f(z - z*) u(<p - y?*) 



(36) 



where 

f(z) = exp 



V2A; 



u{ip) = u(ip,X) 



(1-Xf 



1 + cos ip 



l-2Acosy? + A 2 



(37) 



[note f(z) < 1, u((p) < 1, and /(0) = w(0) = 1]. Here K = V(0, 0) is the maximum 
magnitude of the potential (positive or negative); the center of the patch is at ip = <p*, z = z*, 
Az denotes the axial size of the patch, and A controls its angular size. Indeed, A is related 
to A(p, the angular half-width of the patch [defined as the angle for which u is equal to its 
mean value, u(A<p) = u av ], by the equalities: 



cos Aip = A , A(p = arccos A . 



(38) 



Fourier coefficients of the patch model ( |36i) . ( 1371) are: 



V n (fc) = V„(fc, A, Az) = Vj(k)e-* kz * u n e- m ^ 



1 — A 2 

u n = u n {\) = V2^ A w , n^O; 



f(k) = Azexp 



Uq = U (A) = V27T 



'kAzY 
1 - A 



; (39) 



(40) 



1A ' 1 " ' " " uvv ' "" 2 

For a single patch at each of the two cylinders, the boundary functions G(ip, z) and 
H(ip — 7, z) are given by the formula (13"6"|) as: 

G(<p,z) = V{(p-(pi, X^z-zt, Azx) ; H((p-j,z) = V(y?-p 2 -7, A a , z-z 2 , Az 2 ) . (41) 

The torques corresponding to these distributions are calculated in Appendix B in a closed 
form. Here we study the patch interaction for a particular case when the sizes of the patches 
are identical, Azi = Az 2 = Az, A<pi = Aip 2 = Aip, and V\ = ±V 2 = Vq. 



B. Axial Torque 

As shown in our general analysis, sections III and IV, the axial torque occurs only when 
both cylinders carry non-uniform voltages. For the case of two patches of equal sizes and 
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magnitude, its expression is found by the formula (1B15j) [recall that we have set 7 = 0]: 



7? = VS&z sin 6 (A45) ft sin ( Vl - <p 2 ) exp [-S 2 ] , (42) 

The signs =F in formula fH2l) stand for patches of equal or opposite potential, respectively. 
This torque is proportional to the inverse of the relative gap, d/a. Its dependence on the 
axial patch distance is driven by the Gaussian exponent: the axial torque monotonically 
and rapidly decreases toward zero with increasing values of Z\ — z 2 . The dependence on 
the angular patch distance, < ip\ — ip 2 < tt is shown in figj2]for various values of Aip: the 
torque goes down when the patch angular width gets smaller. 

The axial and azimuthal sizes of the patch play a very distinctive role in the expression 
(l4"2"j) . As Az grows (and the patch thus becomes strip-like, for fixed Aip), the torque mag- 
nitude increases and goes linearly to infinity in Az — > oo. On the other hand, the influence 
of the azimuthal patch width is represented predominantly by the sixth power of the sine of 
Aip: the torque vanishes for belt-like patches, and goes to zero as (Ayj) 6 for Aip — > 0. 



C. Slanting Torques 



1. Uniform and patch potential interaction 

In our case of identical patches expressions (1B3j) for the interaction torque simplify to 

Tl nt = -V^^V V-Azsw 2 {Acp)(zismip 1 Tz 2 smcp 2 ) ; (43) 
Tl nt = V2^^V V~ Azsin 2 (A^)(z lC osv5i =F2 2 cos^ 2 ) • 

Recall that V~ is the difference between the uniform potentials at the boundaries, see 
formula ( flOl) . The minus or plus sign above is taken when the patches have the same or the 
opposite voltages, respectively. Expressions ( 1431) show that this torque is a superposition of 
two contributions each coming from a single patch interacting with the uniform voltage V~ . 
Each of these torques can be expressed as the product of a force acting at the center of the 
patch and the respective arm, z\ or z 2 . It is interesting that these forces are equal to the 
zero order interaction forces obtained in CPEII, section VA. For < Aip < n/2 the torque 
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grows with the azimuthal width of the patch; then, for bigger angular sizes, it decreases 
and goes to zero as the patch becomes annular. Moreover, for Aip — > 0, it goes to zero as 
(Aip) 2 . The torque is again proportional to Az. However, as the patch becomes strip-like, 
the torque goes to zero, by symmetry. In fact, the arm of the force is zero in this latter case. 



2. Patch-patch interaction 



The slanting patch torque expressions ( 1B9[) and ( IBllj) become 
TP = -2v^ € -^V 2 Azsin 2 (Aip) x 



(44) 



M (z\ sin ipi + z 2 sin p 2 ) T ( z i + z 2) M-i (sin ¥>i + sm ¥2) e 



Ty = 2v^ ^-V 2 Azsm 2 (A<p) x 



M (z\ cos <pi + z 2 cos <p 2 ) =F (z\ + z 2 )M \ (cos pi + cos <^ 2 ) e" 



^1 - ^2 
2A^ 



A 



Mi 



1 - A 



1--0 



A) 



1 + A 2 - 2cos(y?i - ip 2 ) 



1 - 2A 2 cos (ip x - ip 2 ) + A 4 



in our case of identical patches. Above we have used a slightly different notation for A^i and 
M as compared to Appendix B, formulas (1B7I) . The first term in the square brackets carries 
the contributions of each single patch independent of their signs. The second term represents 
interaction between the patches. It decreases as the distance between them grows, due to 
the presence of the coefficient Aii and the Gaussian exponent. The decay is faster, the 
smaller the widths Az and Aip are. The plus or minus sign of this term is taken depending 
on whether the patch potentials have the same or opposite signs. Similarly to the previous 
case, the patch torque can be expressed as the sum of the products of a force acting on 
each patch and the arm z\ or z 2 , plus a mutual force between the patches times the arm 
Z\ + z 2 [see figj3] giving torque versus Z\ for different values of z 2 \. The slanting patch torque 
is essentially proportional to the square sine of A(p, and to Az. However, by symmetry, as 
Az — > 00, the torque vanishes, just like the interaction torque above. 



VI. AXIAL PATCH EFFECT TORQUES FOR STEP 

In this section we use the obtained results to evaluate the effect of non-uniform potentials 
on the performance of the instrument that will be used in STEP. The pertinent information 
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about this experiment is found in CPEII, section VIA; for more details, see fl-Q. Here we 
recall the basic design of its core system, the differential accelerometers (DACs). As shown 
in figJU each DAC consists of two test masses (TMs) shaped as coaxial cylindrical shells, and 
of a system of electrodes and magnetic bearings. It essentially constrains the TMs in four 
degrees of freedom, with the translation along their symmetry axis and the spin about it left 
free. The science signal is read by a magnetic SQUID readout system from the differential 
axial motion of the TMs. 

In CPEII we analyzed the patch effect on the axial translation of the TMs. Here we 
concentrate on the axial PE torque causing the spin motion, i. e., relative rotation of the 
cylinders about their common axis. The main reason to examine it is that this motion 
makes PE forces change with the time. Just like in CPEII, we consider a TM and its 
magnetic bearing as a reference case of our pair of cylinders, since the gap between them is 
at least 3 times smaller - and the axial torque thus 3 times larger - than the gap between 
the TM and the electrodes. As was assumed in section |Vj we consider each of the cylinders 
having just one patch, both patches of identical sizes and magnitudes. Even in this case the 
equation of motion proves to be a nonlinear one with the periodic potential such that two 
equilibrium points exist within a single period, one stable and one unstable. We describe the 
general picture of motion pointing out the regimes expected as typical in the real situation 
of the STEP experiment, namely, the regimes of oscillations. We give the expressions for the 
frequency of small oscillations near stable equilibria and the runaway time from the unstable 
ones, and then estimate these quantities under the STEP conditions. An estimate of the 
torque is also provided. Finally, we show that spin motion for any patch voltage distribution 
will be qualitatively the same as in the case of two patches only. 

In compliance with CPEII, we use the following parameters: patch voltage Vq = lOmV, 
TM radius a = 2.3cm, TM to magnetic bearing gap d = 0.3mm, TM length 2L = 0.14cm. 
The moments of inertia for some flight-like TMs vary from 2.34 x 10~ 5 Kgm 2 to 1.40 x 
10~ 3 Kgm 2 for different test masses, but a single TM is fabricated in such a way that the 
principal moments of inertia are all equal to a high accuracy 12|, I x = X y = I z = X. We 
here need, in fact, only X = I z , and we use the smallest value, X = 2.34 x 10~ 5 Kgm 2 . 



14 



A. Spin Motion in the Case of Two Patches 



We take the torque expression (|42|) . and write the motion equation, 

X 7 = T z = =f X* (5) A 2; sin 6 (Ay>) /i( 7 , A) sin 7 , (45) 

for the rotation of the outer cylinder about the symmetry axis; here the dimesionless function 

1 + A 2 

^ = M7, A) = — —-2 > , (46) 

(1 — 2X Z cos 7 + A 4 ) 

and the characteristic value of the torque is 

3 li 

rrt~\ 71 T/ 2 a -z 2 ~ z l — Z 2 { 

Uz) = —e V -e , z = ^^- (47) 

Without any loss of generality, we count here the rotation angle 7 from the position where 
two patches are right one against the other, (^1—^2 = 0. The torque has the minus (plus) sign 
when the patch voltages have the same (opposite) sign. Equation (l4"5j) strongly resembles 
the classical motion equation of the pendulum, with just one additional coefficient /i( 7 , A) 
being a strictly positive non-singular function of 7. Extending this similarity, we note that, 
by the expression ( 1221) for the axial torque through the patch energy W p (j), equation ( 1451) 
has the potential W p {^)/X: it can be equivalently written as 

Z^-« ,48, 

Since the potential is periodic in 7, the complete qualitative picture of motion is well known 



(see, for instance, 



13], Ch.l). It follows from the energy integral of the equation 

^ + W*( 7 ) = E , (49) 
where Eq is the total energy determined by the initial conditions: 

£o = ^ + W p (7o), 7o = 7(to), 7o = 7(*o). 
Potential energy W p {^() is bounded, with the bounds denoted as 

- 00 < W_ = min W p {j) < W + = max W p h) < 00 . (50) 

0<7<27r 0<7<27r 

The minimum potential energy W- = W p { r yJ) corresponds, of course, to a stable equilibrium 
7(i) = 7_ = const, while the maximum one, W + = W p (j + ), is achieved at an unstable 
equilibrium point, 7(2) = 7+ = const. This is enough to qualitatively describe the motion. 
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Indeed, from the energy conservation (1491) it is clear that Eq > W-. If Eq = W-, then the 
system stays at the stable equilibrium, j(t) = 70 = j~, j(t) = 0, t >t . If < £0 < W+, 
then the system can never reach the peak of the potential, the rotation angle is bounded 
at all times, 7_ < j(t) < 'jmax < 7+ ; the motion is finite, which means that the cylinder 
oscillates about the stable equilibrium. If, next, E > W + , then the system always remains 
above all the potential wells, the motion is infinite, the cylinder rotates indefinitely and 
non-uniformly in one direction depending on the sign of the initial velocity. What remains 
is the exceptional case Eq = W + , when, in purely mathematical view, the system stays at 
the equilibrium 7+; however, this rest point is unstable, so in reality any small perturbation 
in this or that direction leads, again, either to the oscillational, or to the rotational motion. 

To make all this even more particular in our case of just two patches, we give the patch en- 
ergy explicitly, as easily found either by the general expression ff2Tl) or, up to an insignificant 
constant, by the direct integration of the torque in the r.h.s. of the equation (H5|) : 

W^) = TUz)^z^{^f) W ( 7 ) + const, ^) = ^ ^ I - 2X^0^ + (51) 

There is no essential difference between the two cases with the opposite signs, so we discuss 
only the case of the minus sign below. 

As seen yet from the equation ( )45l) . in this case we have just two equilibria at each period 
of the potential, the stable one, 7_ = (mod27r), and the unstable 7+ = tt (mod27r). This 
is also clear from the plot of the patch energy given in figj5) the horizontal line through the 
peaks shows the critical energy, W + , that separates the finite motions (oscillations) from 
the infinite ones (rotations). The period of oscillations is the larger, the higher energy Eq 
is, i.e., the larger the oscillation amplitude (the limit of small oscillations, when the period 
is independent of the amplitude, is described in the next section). 

The energy integral (H9l) allows also for the representation of motion in the phase plane 
7, 7. The corresponding plot is given in figj6j closed orbits in it correspond to oscillatory 
(finite) motions (the size of these ovals grows with E ). They are separated from the infinite 
trajectories (rotations) by the so called heteroclinic curves, which go from one unstable 
equilibrium to another nearest to it. The total time of motion along these separatrices from 
one rest point to the other is infinite. 

Note that under the STEP conditions one expects only oscillatory spinning of the test 
mass, rather than its rotation. The reason is that the TM will be caged (fixed) during 
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the satellite launch, and then rather accurately released, practically with no initial velocity, 
which leads to pure oscillations. 



B. Estimates for Small TM Motions Near Its Equilibria in the Case of Two Patches 



Here we describe small motions of the cylinder near its rest points. Accordingly, we 
linearize the equation (j4"5|) by setting 

j(t) =7 T + <5 7 (t), \6j(t)\ < 1 , 

which results in 

h = T <4 *7, 4 = ¥^ sin 6 (Ay?) -&±^) . (52) 



The cylinder thus oscillates about the stable equilibirium position 7_ = with the frequency 



J 2tt 4 M Id sin Ay? ~ sm Ay? 

Accordingly, it rotates exponentially away from the unstable position 7+ = tc with the 
characteristic time 

— —— — * M * ft+«f^)*V > L5 x 10 ,(iwa^ /2 3 (54) 

w+ 7T 3 / 4 V V e a VA7 sin d Ay? VAz sin 3 Ay? 

Numerical estimates (I53p and (I54j) hold for the above set of STEP parameters, with A2 in 
meters. To get the feeling of what the numbers are in reality, let us consider some examples. 
For instance, the longitudinal patch size cannot be larger than the size of the TM, Az < L; 
in the case of the maximum size, the upper bound for the frequency of small oscillations 
becomes: 

< 2.6 x lO-^ 1 + C f ^Hz. 
sin Aip 

This value remains below the STEP signal frequency range 

1.74 x 10- 4 Hz < f < 5.2 x \^Ez ; 

(see CPEII, section VIA) for the patches as small in the azimuthal direction as Aip ~ 1 deg, 
or larger. The frequency tends to infinity when the angular size tends to zero, because the 
potential well becomes infinitely deep. The smallest runaway time in the case Az = L is 



-1/ 



/4 V le a r-r—y/l + cos 2 Ay? _ 5i 2 5 r— a/1 + cos 2 Ay? 



t pe > 5.6 x 10 4 v " ' "" r/ — s 
sin Ay? 
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1 + cos 2 Ay?) 



3/2 



which is much larger than the duration of a single STEP science session, ~ 48 hr, for 
Aip < 10 deg. So if in such case the science session starts with the TM close to an unstable 
equilibrium, then it will be practically stay in this position for its whole duration. 

C. Estimate of the Axial Torque 

By the equation ( )45l) . we got the expression of the axial torque, 

2 1 + A 2 

\T z \=T*(z) Az sin 6 (A(p) sin 7e~ z \ ^(7, A) = — — —3 . (55) 

(1 — 2X 2 cos 7 + A 4 ) 

The ballpark number for the torque is obtained by the maximum of the expression ( 147)) : 

3/2 

^(0) = ^ r eoV 2 ^9.5xlO- 14 iV; 

of the dimension of a torque per unit length. Contributions of the patch widths, and the 
patch azimuthal distribution were not taken into account. A more meaningful estimate 
is obtained by computing the average torque for all the relative positions of the patches. 
Introducing the normalized length Iq = L/2Az, we write, using expression (1551) : 

ITJ = T*Az sin 6 (Aip) / ° dze"^ - / dj ^(7, A) sin 7 , 

l JO 7T JO 

we calculate the two integrals 



l ° _-2 \rK f w 2 
dze z = —erf{l ) , / ^7^(7, A) sin7 = 



lo 2 jo sin Aip (1 + cos 2 A(p) 

and substitute them in the expression, to find 

= ■^T.'^ferfM^^- „ 1.3 x IQ-" (A^er/ik) f ^ JTm 

i/7T L 1 + COS z A<^9 1 + COS^ Ay? 



with Az in meters. For the typical case ^>1, the above expression becomes 

\f z \ « 1.3 x 10- 12 (A^) 2 -g^ Am , (56) 

1 + cos^ Aip 

which is proportional to the square of the axial patch width; for small angular widths the 
mean value goes to zero with the square of Aip, too. 
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D. Spin Motion for a General Patch Distribution 



Here we make an important concluding remark. It is easy to see that the picture of 
the spin motion due to the patch effect torque given in section IVI Al for the case of two 
patches is valid, in fact, for any boundary voltage distributions well. Indeed, in any case 
the electrostatic patch energy is bounded and is a 27r-periodic function of the rotation angle 
7, therefore equations (j48p — (150]) hold, along with the whole following argument about 
the picture of motion based on them. Thus in the most general case the cylinder either 
rotates all the way in the same direction, or oscillates about a stable rest point, depending 
on whether the total energy is above or below the critical value, the global maximum of 
potential energy. The only significant difference is that for a general patch distribution the 
number of equilibria can be larger than two (but always even, with the equal number of 
stable and unstable rest points, since the torque is continuous). 

The increase of the number of equilibrium positions leads to a trend of decreasing the 
amplitude of the typical oscillatory motion. The remark at the end of section IVI Al is also 
valid in the general case: rotational regimes are not anticipated under the STEP set-up 
at all. The experiment conditions are such that any friction or other dissipation should be 
extremely low, so the typical spin oscillations of a TM should be damped but very slowly, as 
compared even with the duration of a single science session. (One should note, however, that 
if some resistive energy losses are present, as it appears to have happened with GP-B 14|, 
then this picture might change significantly). These oscillations will make the axial force 
found in CPEII change periodically with the time. If the basic frequency or some of its 
harmonics is close to the frequency of the STEP science signal, then it might introduce a 
systematic error in the experiment, provided that the PE force and the acceleration due to 
it is large enough. For this reason, as well perhaps for many other ones, careful pre- and 
post-mission calibrations on orbit are recommended, along with the extensive simulations 
before the flight as described in CPEII, section VID. 
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Appendix A: Calculation of the Slanting Torques 



We provide here intermediate results needed for computing of the slanting torques. Par- 
ticular terms that appear under the integrals in the formulas (1251) — (12 8 p are: 



3$p 



COS(f, z- 



smp, z- 



COS<£>, z- 



sin cp, 



COS (f, 



smp ; 



(Al) 



dp dp dp dp dz dz 

we need to determine the Fourier coefficients of these functions evaluated at the boundary 
p = b. The radial and the angular derivatives of the potential involved here are: 



dp 



P =b 



d$p 



dp) 



1 x 

d 



2tt 



dk J2 



z(kz+nip) 



G n (k) - H n {k)e-^ 



dk ^ < 



i(kz+nip) 



in H n (k)e 



-4717 



(A2) 



(A3) 



they come from the formulas (fT9|) and ( !3T|) . respectively. The z derivative is computed from 
the second of the boundary conditions ( TTTj) : 

d$p 



dz 



^i(kz+rup) 



p=b 27T J -°° n=-oo 



kH n (k)e- 1 ^ 



(A4) 



For any function u(p, z) such that zu(p, z) is squarely integrable, by the definition of the 
Fourier transform ([7]), the following equalities hold: 



2 roo 00 



du n (k) _ 
dk ' 



(A5) 



1 />00 °° 

u{<p,z)co8<p = — ]T e« kz+ ^ [u n ^(k)+u n+1 (k)} ; 
4 7 rJ-oo n= _ 00 

• CO 

u(p,z)smp = - — / Yl e i{ - kz+nip) [u n -i(k) - u n+ i(k)] . 

Using formulas (IA5I) and (IA2I) with the radial derivative playing the role of u(<p, z), the first 
two functions (1A1I) are represented in the following way: 



z— — cos cp 
dp 



P =b 



And 



dk J2 



i(kz+nip) 



d_ r 
dk 



(G n ^(k) - H n ^(k)e-^-^) + (A6) 
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z— — sm 09 
dp 



And 



OO Q 

dk ei{kz+nip) ^ 



(G n+1 (k) - H n+1 (k)e-« n+1 ^) 
G n . x (k) - H n ^{k)e-< n -^) - (A7) 
(G n+1 (k) - H n+l (k)e-^ n+1 ^ 



By the same token, using formula (1A3I) instead of (1A2j) . for the two terms proportional 
to the potential derivative in 09 we obtain: 



— COS 09 

dip 



z— — sin 09 

dip 



7Z dk H 



p=b 



47T 



I l-oo 
An 



d_ 

dk 



'n-l)H n ^{k)e- i{ - n -^+ (A 



/OO x 
dk y: 
-OO , . 



x d 

i(fc2+ny3) 

dk 



n + l)H n+1 (k)e- i( - n+1 * 
n+l)H n+1 (k)e- i( - n+1 * 



(A9) 



Note that expressions (1A6j) — (1A9j) here are valid under the additional conditions 



5 n=— oi 

which are equivalent to: 
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(A10) 

To determine the remaining two terms in (lAlj) containing the derivative of the potential 
with respect to z, we just need the last two formulas from (1A5|) and the expression (1A4j) in 
place of u(cp, z). This results in: 
d$P 



dz 
d$P 
dz 



COS 09 



sm 09 



7 yoo ^ r -, 

/ dk y e i{kz+nv) k (fc)e- i(n - 1)7 + ^ n+ i(A;)e- i(ri+lh ; (All) 

47T J-oo 



1 

47T j-x 



ra=— oo 
oo 



dk Y e 4(fc2+nv) fc H^ky-^- 1 ^ -H n+ i {k)e~ i{n+iy< . (A12) 



Formulas f ]A6j) — fl Al 2|) provide the expressions and the needed conditions (lAlOj) for cal- 
culating the integrals in the formulas (125 p — (I28p for the slanting torque. 



Appendix B: Calculation of the Torque for a Single Patch at each of the Cylinders 

To get the torque we first need to find the Fourier coefficients of the boundary distributions 
G(ip, z) and H{tp — 7, z) from the formulas (jUJ) and their derivative with respect to k. The 
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former can be represented by the Fourier coefficients (|39|) . (|40j) as: 



G n (k) = V n (k, Ai, A**) ; if n (fc) = V„(fc, A 2) A^ 2 ) 



(Bl) 



The latter are the derivatives in k of these expressions, and they are given by the general 
formula: 

2" 



dVj(k) 
dk 



-VjAz [k(AzY + izj) exp 



/ kAzj 



MA;) e- i(fc *' +n ^ ; j = 1,2 , 



(B2) 



with V*(Jfe) = G n (Jfe) and V^ife) = H n (k). 

The torque due to the interaction between the patches and uniform potential difference 
is a linear function of the derivatives in k of G n (k) and H n (k) calculated at k = and n = 1 . 
These are determined by the formulas (IB2I) as: 



dG x {k) 



dk 



l-\\ 



-vpi . 
1 



dGx{k) 



k=0 



dk 



2TfV 2 Az 2 (iz 2 



k=0 



where we used equality ( 140)) to express the coefficient u n . In order to calculate the torque 
we need just to substitute the above expressions in the formulas lj32l and (|33l) . This leads 
to the following result: 



T^ = -V2^ 6 -^V- 
d? 



ViAzi (sin Aipi) 2 z\ sin<£>i — V 2 Az 2 (sin A<£> 2 ) 2 z 2 sin ((p 2 + 7) 



V\Az\ (sin A(/?x) 2 z\ cos^x — V 2 Az 2 (sin A(p 2 ) 2 z 2 cos (<p 2 + 7) 



(B3) 



Note that here we used the first of the relations (1551) to replace A with the more meaningful 
parameter Atp . 

The expressions for the slanting torque due to the interaction between the patches are 
more cumbersome to find, since one needs to calculate the sum over n and the integral over k 
of the product of the Fourier coefficients of the boundary distributions and their derivatives 
in k. For the x component, we combine formula (134j) with (1B2I) for the derivatives, and ( TB1I) 
for the boundary functions, to obtain the expression: 



= 2vr ^ r dk m -VfAzlN^) \kAz\ + izA exp \-k 2 Azl 

d J — 00 L J L 



-lipi 



+ 



(B4) 



V 1 V 2 Az 1 Az 2 M 1 



kAzj + iz 2 ) e lk{zi - Z2) + (kAzl + iz x j < 



-ik{z\—Z2) 



exp 



Az\ + Az\ 
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V 2 2 Az 2 2 N(X 2 ) kAz 2 2 +iz 2 



exp 



-k 2 Az 2 



-Uf2 



for the matter of space we set 7 = 0, without any loss of generality, and denoted [compare 
with CPEII, Appendix A]: 

00 



M 1 = M 1 (X uV > 1 ,X 2 , V > 2 ) = ^- £ u n (X 1 )e^u n+1 {X 2 )e-< n+1 ^ ; 



(B5) 



00 

N{\) = — u n (X)u n+l (X) = M 1 (0,X;0,X; 



(B6) 



The values of these coefficients are obtained by summing up geometrical progressions [see 
coefficient u n in formula (1401)]: 

e -*((Pi-!«) _ A X A 2 



Mi 



1-Ai 1-A 2 



1 + A 2 



l + y(l + Ai) 



l + y(l+A 2 , 



+ 



(B7) 



.0 = 1- 2AiA 2 cos(y?i - <£ 2 ) + (AiA 2 )' 



N(X) = 1-^(2 - A) . 



(B8) 



All we need now to get is the two integrals in the formula (IB4j) . which are well known: 



d/c exp 



k 2 (Az 2 + a^ 2 ; 



'2k 



Az 2 + Az 



exp 



2 (Az 2 + Azf) 



(ifcexp 



fc 2 (Az 2 + A^ 2 ; 



ke ±lk ^- z ^ = ±iV2^ — ^^exp 

(Az 2 + Az 2 2 f 2 



,zi - z 2 



2 (Azj + Azl) 



In the case z\ = z 2 and Az% = Az 2 = Az they becomes: 



dk exp 



-A; 2 Az 2 



\px I Az ; / dkexp 



-k 2 Az 2 



k = . 



With these results we are now able to rewrite formula ( 1B4I) for the torque as: 
= -2tt 3 / 2 ^jyfc A3xiV(Ai) sin^x + V r 2 2 z 2 Az 2 iV(A 2 ) sin^ 2 - 

l/ x l/ 2 Z~5 (Mi) (Az 2 - Az 2 ) + ( Zl + z 2 )] exp [-z 2 ] } , 



(B9) 



with the new notations 



/ = V2 



Az t Az 2 
Azl + Az 2 



Zl - z 2 



2 {Azl + &zl 



(BIO) 
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The calculation of the y component of the patch torque does not present any additional 
difficulties. In the same way as we derived formula flB9j) . one can find for T^: 



= 2vr 3 / 2 ^ ju^A^Ax) cos^ + V 2 2 z 2 Az 2 N(X 2 ) cos<^ 2 
V X V 2 m (Mi) [y/2z (Azl - Az\) + ( Zl + z 2 



(Bll) 



exp [— z 2 } 



[again, the notations (1B7|) and (IBlOj) are used]. 

To obtain the closed form representation of the axial torque we combine the formula ( )3T|) 



with the expressions fIBlj) . The integral in k there is found above, so the result is: 



7J = 4vr 3 / 2 ^U X U 2 / 9 (M») exp [-5 2 ] , 



(B12) 



with the coefficient M 3 given by 



1 00 

M 3 = M 3 (\ 1 ,<p 1 ,\ 2 ,<p 2 ) = — £ n Mn (A 1 )e- i ^ 1 Mn (A 2 )e^ 2 

Z7T ' 



(B13) 



Formula ( I4TJ]) for u„ leads to an explicit sum of this series reduced to the derivative of a 
geometric progression: 



M, = -i 



;i-A?)(l-A 2 



a - a i a 2) 



sin (<y?i - v? 2 ) 



(B14) 



Using this we write formula (1B12I) in the final explicit form: 



3/2 

Tl = ^ Vl V 2 I (l - Xfj (l - A 2 ) sin fa - V2 ) 



(B15) 
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FIG. 1: Geometry of the problem and coordinate systems 
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FIG. 3: Slanting torque vs. the axial distance between the patches for Z2 = , Az , 5Az. 




FIG. 4: Cross section of the STEP differential accelerometer (all dimensions in meters) 
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FIG. 5: Periodic potential 
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6: Phase plane of spin motion 
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